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The Maxwell's equations are solved when it has an inhomogeneous terms as a source. The solution is very 
general in a sense that it handles arbitrary current source and anisotropic media. The calculation is carried 
out in the k-domain after Fourier transform, and its results are confirmed with the propagators of the new 
coupling-free wave equations derived from the Maxwell equations for anisotropic environments. . 
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The set of Maxwell's equations is a classical repre- 
sentation of electromagnetic field theory. It is still a 
center of technological applications as well as theoretical 
research in many areas of physical science. The physical 
prediction of wave propagation is made through the 
fields calculation and its basis is the Maxwell's and 
wave equations. The task obtaining the solution of the 
Maxwell's equations is important as much as under- 
standing the the equations themselves. 

The Maxwell equations are the first order, linear, 
simultaneous partial differential equations (PDE) for 
vector fields E and H, but their solution is not easily 
obtained. The main difficulty lies in the fact that the 
unknown variables E and H are coupled to each other 
in the equations. The most frequently cited method of 
solving the Maxwell's equations is to take a curl opera- 
tor on both sides of Maxwell's 3rd(Faraday's law) and 
4th(extended Ampere's law) equations. This operation 
is specially powerful for an isotropic medium, since it 
gives a decoupled version of wave equations with source 
terms. These wave equations are not only the equations 
for waves, but also solvable forms of the Maxwell's 
equations. The wave equations of the isotropic medium 
have been solved conventionally by Green's function 
technique. 

The solution of the wave equations has a long history. 
The practical advancements have been achieved in 
1950s by Schwinger and Marcuvitz, who solved first the 
inhomogeneous wave equations using Green's functions 1 . 
Soon after, Many other scholars, for example, F. V. 
Bunking C. T. Tai^, and L. B. Felsei4 developed 
methods using the so-called dyadic Green's functions, 
which have been used to this time. 

The tensor(or dyadic) Green's functions have been 
widely studied specially last decade, as the number of 
electronic equipment operating in an anisotropic envi- 
ronment increased. However, the general solution of 
Maxwell's equations or wave equations in an anisotropic 



medium have not been known because of coupling prob- 
lem between E and H. The coupling problem of the equa- 
tions for the anisotropic medium is more complicate than 
the isotropic case, since the permittivity e and perme- 
ability p are matrices. The general wave equations for 
an anisotropic medium become 
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It is noted that the first term of Eq.(l) was not expanded 
to V(V • E) — V 2 E, as the V • E cannot be written as 
iT 1 p when applying the Gauss's law. The correct form 
of the first Maxwell's equation is V • eE = p, instead of 
V • E = e p, since e is a matrix. The e cannot move 
out from the inside of divergence. It is also important to 
point out that the e and p in the curl operator of Eq.(l) 
and (2) cannot come to the front of the curl operator, 
as a constant number does. Thus, the equations can be 
solved only within a restricted range. In an arbitrary 
general case, the equations are practically regarded as 
insolvable without using numerical simulation as far as 
we follow usual conventional ways. 

This paper attempts general solutions of the fields E 
and H in 2 different ways. The first method is to solve 
the fields directly from the Maxwell's equations, and the 
other one is to apply the Green's function technics on 
the new wave equations, which are derived in a different 
way to fit into an anisotropic environment. The two 
solutions, which are obtained respectively in different 
ways, turn out to be same. 

It is useful to write down the Maxwell's 3rd and 4th 
equations in an unifying description method as following, 
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By using harmonic time dependency, the above 
equation becomes 
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E, H, and J in this equation are now the function of 
spatial coordinates. The equation is separated in 3 parts 
- operator, unknown E and H, and inhomogeneous term 
as a source. The operator part in the bracket acts on 



fields 



which constitutes a linear simultaneous 



PDE. The Fourier transformation has an advantage 
in transforming a linear differential equation into a 
corresponding algebraic equation. The curl operator 
is replaced by (— zkx), but this vector notation is 
not convenient in conjunction with matrix structure. If 
we use the following matrix k representing wave vector k 5 



(5) 



then k x A = k ■ A for any arbitrary vector A. There- 
fore the Eq.(4) becomes 
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Here, k is defined in rectangular coordinates. This im- 
plies that k is coordinates dependent. One can find it 
easily in spherical or cylindrical coordinates. The quanti- 
ties E, H, and J in the equation are Fourier transformed 
functions of k. The first 2x2 matrix is actually a coeffi- 
cient matrix of the simultaneous equation, but it should 
be taken care of the fact that each component of matrix 
is 3 x 3 matrix. The 3rd and 4th Maxwell's equations be- 
come simultaneous equations. The solution is calculated 
by diagonalization of the coefficient matrix. 
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As E and H are functions of k, the original spatial 
fields E(r, r') and H(r,r') are obtained by an inverse 
Fourier transformation. The spatial points r, r' are each 
point of field and source. Now, the solutions are derived 
directly from the Maxwell's equations, not from the 
wave equations. Moreover, the application range of the 
solution is extended to an arbitrary anisotropic media. 

As seen above, the simple algebraic manipulation of 
Maxwell's equations leads to the solutions. If one takes 



look at the Eq.(7) carefully, the right term [ ) 

plays a role of source term, and the components of the 
first matrix in the left side can be interpreted as Fourier- 
transformed Green's functions or 3 dimensional propa- 
gators of the fields. Usually, the anisotropic propagator 
of the electromagnetic fields are available only in special 
cases, but the solution (8) and (9) provides the general 
propagators as follows. 



g E (k) = (k- ^k + ^e)- 1 
gH^^Ck-e-^k + Lo 2 ^)- 1 
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The fact that we find the general propagators implies 
we can calculate fields only by integration. Thus, the 
solving procedure of the electromagnetic fields can be 
automated. 

If these solutions, directly obtained from the Maxwell's 
equations, are right, the wave equations should also 
give the same results. However, e and n are replaced 
with matrices in the anisotropic medium, and the wave 
equations cannot be decoupled in general. The problem 
is essentially due to the fact neither e nor fi can be 
extracted from the vector differential operators. Unlike 
an isotropic case, it is impossible to apply Maxwell's 
equations on V x eE and V x /iH. The terms e • E and 
H ■ H are merely a linear combination of the original E 
and H, but they act as another unknown independent 
variables, when contained inside a vector differential 
operator. With Eq.(l) and (2), the decoupled version 
of wave equations is not possible. There are some 
examples of solvable models even in the anisotropic 
medium, when the anisotropic characteristic comes from 
either e or /i 6 . However, the general problems, when 
anisotropy occurs due to both e and /z, can not avoid 
the coupling problem. Whether they are solvable or 
not depends on whether they become decoupled or not. 
The solution of the general problems doesn't seem to 
be feasible as far as the ordinary wave equations are used. 

At this stage, the goal is obviously to find the de- 
coupled version of wave equations capable of solving the 
fields E and H. A better approach is to go back to the 
original Maxwell's equations and derive new equations 
instead of relying on the original wave equations. The 
magnetic field H can be written from the 3rd Maxwell's 
equation as 
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(V x E) dt. 
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Inserting the above equations into the 4th Maxwell's 
equation and differentiating with respect to time gives 



V x (/x" 1 • (V x E)) + e 
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This equation is a new type of a wave equation for the 
field E. It is decoupled and doesn't have H. A new wave 
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equation for H can be obtained in a similar way. The 
result reads 



V x e _1 (V x H) 
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Now, these two equations Eq.(13) and (14) constitute 
new wave equations^ instead of the original wave equa- 
tions for the isotropic case. These new equations are de- 
coupled even in any arbitrary anisotropic environment, 
and consequently solvable. Moreover, they recover their 
original wave equation forms when there is no anisotropy. 
The next step is to check whether or not the propagators 
of the new and the original wave equations are identical 
to each other. The time dependency is again assumed to 
be harmonic. Then the propagators of Eq.(13) and (14) 
are obtained from the spatial Green function equations 
replacing source terms with the S functions. The Green's 
function Ge(^,t') and G^(r, r') are 

V x (/I" 1 ■ (V x G E )) - u?t G E = S(r - r'), (15) 

V x (e^ 1 • (V x G H )) - 0J 2 fi G H = S(r - r'). (16) 

One of the advantage of Fourier transformation is in 
switching a linear differential equation to an algebraic 
equation. The propagator is the Fourier Transform of 
spatial Green's functions. 

G e (t,t')= J ff£ (k)e-* k '( r - r Vr (17) 

G H (r,r') = J 5ff (k)e-' k -( r - r Vr (18) 

The k-space algebraic equations for propagators gg (k) 
and gn (k) are obtained by inserting above equations into 
Eq.(13) and (14). 



g H (k) = -(jee^k + uPe)- 



(19) 
(20) 



Comparing with the results obtained in the previous 
paragraph, there exists (-1) times difference overall in 
both gE and gn- These differences come about due to 
(-1) times differences in the source terms. The both 
results are exactly same. 

The two equation sets, Maxwell's equations and wave 
equations, which look different in the r space, turn out to 
be same in the k space. The fact that two equations have 
the same solutions has an mathematical implication. 
The solution of Maxwell's equation has been obtained 
through the diagonalization of the matrix, but the other 
solution of the wave equation has been treated with 
the Green's function technique. The coincidence of two 
solutions indirectly shows supports through the mutual 



verification of both solutions. 

As for the justification or the proof of the solutions, a 
few things can be mentioned. Firstly, the new wave equa- 
tions become the original wave equations, when there is 
no anisotropy. It implies that the new equations are the 
extension of the original wave equation to the anisotropic 
environment. Secondly there are a few number of re- 
searches, whose results are consistent with that of this 
article. For example, Cottis, Vazouras and Spyrou^ cal- 
culated the Fourier-transformed dyadic Green's functions 
in an electric anisotropy, and their result reads 
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where e r , and /i r are the relative matrix coefficients as 
usual. The result of the new equation for an identical 
case is given below according to Eq.(19): 



gE NEW (\t) 



-1 



(22) 



Given that there is no anisotropy in fi, the identity 
matrix is substituted in place of (i^ ■ The same result 
is obtained using the relationship k ■ k — k ® k — k 2 I. 
In addition, calculation of magnetic anisotropy and 
conditions of refractive index imposed on plain wave 
also show identical results between two systems. This 
coincidence is natural consequences, because they are 
based on the same mathematical ground. 

In conclusion, we have shown the general solution of 
Maxwell's equations with source terms in 2 different 
ways. The solution of the first attempts was calculated 
directly from the Maxwell's equations using matrix diag- 
onalization of operator part after Fourier transformation. 
The second method is to solve the newly derived cou- 
pling free wave equations with the tensor Green's func- 
tions. The two solutions, from matrix diagonalization 
and Green's functions, turn out to be identical, and are 
coincident with previous researches. These solutions are 
expected useful in many practical field calculations in- 
cluding inverse problems. 
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